Fluctuation theorems in inhomogenious media under coarse graining 
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We compare the fluctuation relations for work and entropy in underdamped and overdamped 
systems, when the friction coefficient of the medium is space-dependent. We find that these relations 
remain unaffected in both cases. However, for the overdamped system, the analysis is more involved, 
and a blind application of normal rules of calculus would lead to inconsistent results. 
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I. INTRODUCTION 

The last couple of decades have observed a steadily 
growing interest in the field of systems at mesoscopic 
scales, thanks to the growing understanding of machines 
and engines with smaller dimensions. This has led to 
the area of stochastic thermodynamics which provides a 
framework for extending notions of classical thermody- 
namics to small systems wherein concepts of work, heat, 
and entropy are extended to the level of individual trajec- 
tories during nonequlibrium processes (ensembles). Re- 
search in this area has given birth to a group of exact 
and powerful theorems that dictate the behavior of such 
systems. They are commonly referred to as the fluctua- 
tion theorems (FTs) |ll4l6l| , and these theorems are valid 
even far from equilibrium, a feat that is beyond the scope 
of the well-established linear response theory. The theo- 
rems provide stringent restrictions on the probabilities of 
phase space trajectories in which second law is transiently 
"violated". They show that at the level where fluctua- 
tions are comparable to the relevant energy exchanges of 
the system, one needs to replace the associated quanti- 
ties in the statement of the second law by their averages: 
(W) > AF or {Astot) > [1 [H, [ll] . Here the angular 
brackets represent the ensemble average. Thus, they in 
essence uphold the second law, even at the mesoscopic 
level, however, for the average properties. 

The Crooks Fluctuation theorem (CFT) for heat states 
that the ratio of the probabilities of forward trajectory 
and the correspondin g re verse trajectory for given initial 
states is given by [l^ [la | 



P[X\xo] 

PlXlXr] 



„/3Q 



(1) 



Here, X is the short form of the phase space trajectory 
along the forward process Xq,Xi, ...,Xt generated by the 
protocol A(i). Xi represents the phase space point at time 
ti. X is the corresponding reverse trajectory generated 
by the time reversed protocol Ar-t, where r is the time 
of observation, xo is a given initial state of the forward 
process. The reverse process begins from the state Xr, 



which is the time-reversal of the final state Xj- of the 
forward process. 

Using CFT, several other theorems like the Jarzynski 
equality and entropy production FT, can be easily de- 
rived [il[ll]. 

In this paper, we study the validity of these FTs in the 
presence of coarse-graining . We flnd that a prominent 
difference in the analysis is observed between the over- 
damped (coarse-grained) and the underdamped systems, 
when the friction coefficient is space-dependent '17H20| . 
We find that although the derivation is simple for the un- 
derdamped case, the overdamped case will lead to incon- 
sistent results, if normal rules of calculus are blindly ap- 
plied to it. It should be noted that space-dependent fric- 
tion does not alter the equilibrium state vis-a-vis space- 
independent friction. However, Langevin dynamics of 
the system gets modified especially for the overdamped 
case. There are several physical systems wherein friction 
is space-dependent (see [20] and the references therein). 



II. CROOKS THEOREM IN PRESENCE OF 
SPACE-DEPENDENT FRICTION 

In the presence of space-dependent friction 7(2:), the 
equation of motion of the underdamped system of mass 
m moving in a time-dependent potential U{x,t) is given 
by 



mv = -j{x)v - U'{x,t) + yj2'j{x)T£,{t). 
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(2) 



This can be derived microscopically [18|, [19| . Note that 
the above equation contains multiplicative noise term. 
Here, T is the temperature of the bath, while ^(i) is the 
delta-correlated Caussian noise with zero mean: (^(i)) = 
0; (^(i)C(i')) = Ki - t')- The overhead dot denotes time- 
derivative, whereas prime represents space derivative. 
Eq. ([2]) has been derived microscopically invoking system 
and bath coupling [l^, [l^. It is shown that the high 
damping limit of eq.© is not equivalent to ignoring only 
inertial term [13-120] • The detailed treatment leads to an 
extra term that is crucial for system to reach equilibrium 
state in absence of time-dependent perturbations (see eq. 
[TO] below). Let us now check the validity of CFT in both 
the underdamped and overdamped cases. 



A. Underdamped case 

At first we want to calculate the ratio of path prob- 
abilities between forward and reverse process. In a 
given process, let the evolution of the system in phase 
space be denoted by the phase space trajectory X{t) = 
{xojXi, • • • ,Xr}- Here, Xk represents the phase point at 
time t = tk- In general, the phase point includes both 
the position and the velocity coordinates of the system. 
In the overdamped case, however, it would consist of the 
position coordinate only. Now, a given path X{t), for a 
given initial point xq, would be fully determined if the 
sequence of noise terms for the entire time of observa- 
tion is available: ^ = {^OiCi: ' ' ' iCt-i}- The probability 
distribution of S,k is given by 



P{^k) oc e-^>^\ 



(3) 



Therefore, the probability of obtaining the sequence ^ 

will be fia, m 



P[^{t)] oc exp 



1 



C^{t)dt 



(4) 



Now, from the probability P[^{t)] of the path ^(t) in 
noise space, we can obtain the probability P[X(i)|xo]. 
These two probability functionals are related by the Ja- 



cobi 



I ail 



Thus, we can as well write [1/ 



P[X{t)\xo] ocexp 



f{t)dt 



(5) 



where the proportionality constant is different from that 
in eq. ((4]). In eq. ([5]), we then substitute the expression 
for ^(i) from the Langevin equation (Eq.©): 



P[X{t)\xo] ocexp 



dt 



{■mv + U'{x,t) +-/{x)vy 



j{x)T 



(6) 

For the reverse process, v — >■ —v, but the Jacobian is 
same. The ratio of probability of the forward to the re- 
verse path can be readily shown [I^, [22| 

PlX{t)\xo] 



P[X{t)\Xr] 

_ exp [- Jg dt{mv + U'{x, t) + 'y{x)vf /A'y{x)T\ 
^ exp [- j^^ dt{mi + U'{x, t) - ')(x)vY /'^')(x)T] 



— exp 
= exp 



dt 



4m'y{x)vv + 4U'{x, t)^{x)v 



4j{x)T 
P I dt{mvv + U'{x,t)v) 



(7) 



where Q is the heat dissipated by the system into the 
bath, defined as 



{-r{x)v - ^2j{x)TC{t)}v dt 
{mv + U'{x, t)}v dt. 



(8) 



This definition follow from the stochastic energetics de- 
veloped by Sekimoto [23i, [24] from the definition of First 
Law using Langevin dynamics . Eq.([7|) is the celebrated 
CFT, from which several FT follow. 



B. Integral and detailed fluctuation theorems 



We have. 



P[X(t)\xo] 

P[X{t)\ir] 



aPQ 



(9) 



where Q is the heat dissipated, as obtained from the first 
law. Multiplying by the ratio of the initial equilibrium 
distributions, for forward and reverse processes, namely 
by Po{xo)/pi{xt), we get [ll| 



P[X{t)\xo]po{xo) 

P[X{t)\S:r]PliXr) P[X] 



W^e^Q 



-/3-Bo 



Z{Xo) 



Z[K) 



8{Q+AE~AF) 



„/3(W-AF) 



(10) 



We have used the expression for equilibrium initial dis- 
tribution Pq{xo) — '^.^ ^ and Pi{xr) = ^zUT)' -Here, 
AE = Er — Eq, and we have made use of the relation 
Z — e^^^, between the partition function and the free 
energy. Z{Xq) and Z{Xt) are the partition functions cor- 
responding to the protocol values at the initial time and 
the final time, respectively. In the final step, the first law 
for the work done on the system, W — Q + AE, has been 
invoked. The above relation can be readily converted to 
the Crooks work theorem [i6| , given by 



P{W) 
P{-W) 



„/3(W-A_F) 



(11) 



Where P{W) is the probability of W work done on the 
system in the forward process, while P{—W) is the prob- 
ability oi —W work done on the system in the reverse 
process. By cross-multiplication and integration over W, 
we get the Jarzynski equality 9]: 



( 



^-0W^ ^ ^-0AF^ 



(12) 



If the initial distributions for the forward and reverse 
processes are not equilibrium ones and pi{xt) is the so- 
lution of the Fokker-Planck Equation at the final time r 
of the forward process, we get, instead of eq. pU]). the 
relation [11, [H 



P[X] 
P[X] 



pl3Q+\n{po(xo)/pi{x^)) 



a^Sto 



(13) 



We then arrive at the relations for change of total en- 
tropy Astot which is nothing but sum of change of system 
entropy Assys = \n{pQ{xo)/pi{Xr)) (in the units of Boltz- 
mann constant ks) and entropy production in the bath 

SB = PQ- 



Astot = ln{poixo)/pi{xr)) + PQ- 



(14) 



Tfcl.O 




0.05 0.1 0.15 0.2 



FIG. 1: Transient work distribution for underdamped case 
with space-independent friction 
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FIG. 2: Transient work distribution for underdamped case 
with space dependent friction 



From Eg. (1131) integral fluctuation theorem follows, which 
hold for all times, namely. 



(e-^^-0 = 1- 



(15) 



From the integral forms of the fluctuation theorems, 
given by eqs. (IT^ and (ITSI) . using Jensen's inequality we 
obtain easily second laws jOi, [ijj 



{W) > AF; 
(Astot) > 0. 



(16) 

(17) 



Thus, in the underdamped limit second law retains same 
form for a system in presence of space-dependent fric- 
tion. This completes our treatment for some FTs in the 
underdamped case for a particle moving in space depen- 
dent friction. 



Above exact FTs do not give any information 
about probability distribution of work {P{W)), entropy 
P{Astot) etc. These distributions depend crucially on 
the specific problem being investigated. 

Here, we study these distributions for the case of driven 
particle in harmonic trap. Apart from verifying FTs we 
also see how the space dependent friction modifies the 
distribution of {P{W)), and P{Astot) as compared to 
the particle moving in a space independent frictional co- 
efficient 7 (which is the space average of 7(a;)). The un- 
derdamped Langevin equation is given by 



mv = -7(x)u -kx + A sin(wt) -I- y/2-/{x)T^{t). (18) 

Asm{ujt) is driven sinusoidal force of frequency uj and 
amplitude A. For this model analytical solution can be 
obtained for space independent case only for both over- 
damped and underdamped case [31 , 32] . 

For simplicity in our study, we restrict ourselves to two 
cases of space dependent friction (i) 7(x) — j — constant 
(ii) j{x) = 7 + ctanh(aa;) 

In fig.© we have plotted the transient work distri- 
bution obtained after driving a system for one-fourth of 
a cycle for forward {P{W)) and corresponding reverse 
{P{—W)) protocol. Initially the system is equilibrated 
at appropriate initial values of protocol for forward and 
reverse process. In all our simulations, we have used 
the Heun's method of numerical integration |25|, and 
have generated '^ 10^ realizations. Implementing the 
Heun's method tantamounts to using the Stratonovich 
discretization scheme [2a]- Henceforth we have used all 
the quantities in dimensionless form and taken k=l, m=l 
and 7 = 1. For case (i), both distributions are Gaussian 
nature, and they cross each other at AF = —0.044, which 
is the free energy difference over one-fourth cycle. This is 
obtained numerically from (e"'''^) = e~^^^, while the- 
oretically we have AF = —0.045. This is well within our 
numerical accuracy. 

In fig© we have plotted the same for space dependent 
friction 7(x) — 1.0+0.9 tanh(20a;). Here the distributions 
are non-Gaussian but the crossing point is same as in 
space independent case. This is because the equilibrium 
distribution remain same in both cases. 

In fig. ([3]) we have plotted distribution of total entropy 
production of a driven Brownian particle confined in a 
harmonic trap for one- fourth cycle. Here the underlying 
dynamics is underdamped and we find that the distri- 
bution is Gaussian for space independent case, while it 
is non-Gaussian for space dependent case. If we take 
particle to be initially equilibrated at different tempera- 
ture T = 0.05 and then connected instantaneously to the 
given bath of temperature T — 0.1, and driven by same 
external force (i.e, for athermal case), we find that the 
distribution of total entropy production is non-Gaussian 
even for space independent case. This is consistent with 
the results in ^. Numerically we find (e^^"'"') = 1.002 
which is well within our numerical accuracy. In all these 
distributions, we find that there is a finite weight for re- 
alizations having W < AF and Astot < 0, although the 



7(x)^l+0.9tanh(20x) - 
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FIG. 3: Distribution for total entropy production for under- 
damped dynamics 



mean values follow the second law inequalities. These re- 
alizations are called transient Second Law violating tra- 
jectories. This finite weight is necessary to satisfy the 
fluctuation theorems j28l |. 

After establishing the well known FTs in the under- 
damped case, we turn our attention to the overdamped 
dynamics of the particle, in a space-dependent fric- 
tional medium. Going to the overdamped regime implies 
coarse-graining. Instead of evolution in full phase space 
(coordinates and momenta), we restrict the evolution of 
the system to the position space only. This is equiva- 
lent to ignoring the information contained in the velocity 
variables. 



C. Overdamped case 

The treatment of overdamped case is more subtle and 
need to follow proper methodology. The overdamped 
equation is modified depending on the discretization pro- 
cess and it can be written as 



-r{x)u'{x,t) + {i 



a)9{x)9'{x)+g{x)^{t). 



For detail we refer to [20| . Such ambiguity of discretiza- 
tion process does not arises in the underdamped case as 
discussed in detail in [23|- Here, g{x) — ^/2T^{x) = 
^J2T l^(x). a G [0, 1]. a = for Ito convention, while 
a — Xjl and a — \ for Stratonovich and and isothermal 
conventions, respectively. In [201, i^ ^^^ been shown that 
for all values of a, the same equilibrium distribution is 
obtained for a given value of the pro tocol. Now we follow 
the treatment closely as given in [20| . For a given eq. (fT9)) 
the path probability for a single trajectory in position 
space is given by 



P{X[i)\x^\^e 



^s\x\ 



(20) 



where 



5[X]=^ dt[^\x-!{x,t)^agg'f 



+ af{x,t) 



(21) 
Using f{x,t) = -U'{x,t)T{x) + (1 - a)g{x)g'{x), we get 

S[X]= j^ dt(^[x-U'T+{2a-\)ggf 

+a[-U"T - U'V + (1 - a){gg" + g'^)]) . (22) 

For reverse path, (see eq. (22) of [30|,[33|) one has to 
replace x — >■ — i, and a — )■ 1 — a. Thus the action for 
reverse path is given by. 



S[X] 



+ il-a),f'{x,t)). 



(23) 



Once again, substituting f{x,t) = —U'{x,t)T{x) 
ag{x)g'{x), we get 



1 



S[X] = I dt[^^[-x- U'T - (2a - l)gg'] 

„'2 



\2g' 

+ (1 - a)[-U"r - U'V + aigg" + g'')]) . (24) 



Using eq. ([22|) and eq.(|24|) the ratio of forward and reverse 
trajectories will then become 



P[X\xo] 

P[X\Xr] 



^^SIX]-S[X]. 



(25) 



(19) 



where 



S[X] - S[X] 



dt 



2g^ 



■[-X - U'T - (2a - l)ggf + (1 - a)[-U"T - U'V + a{gg" + g'^)] 



-4iC/T - 4(2a - l)U'gg'T 



= / d< 



= / dt 



dt\ 



252 
-'ixU'T ~ 4(2a - l)C/'rrT 



4Tr 



(2a-l)(C/'T + C/T') 



+ (2a-i)(c/"r + c/'r') 



<it[ 



-/3iC/' - (2a - l)C/'r' + (2a - \){U"T + C/T')] 
-l3xU' + (2a - l)[/"r] , 



(26) 



where we have used g^ = 2TV, and gg' — TV. From the 
generahzed chain rule in presence of muhiphcative noise 
[sdl we have 



^^(^'*)=f + t/'x+^(l-2a)t/" 



dt 



(27) 



It may be noted that — ^^ is total rate of change in- 
ternal energy, whereas ^dt is the differential work done 



on the system. Thus, 

S[X] - S[X] = -/3 




/3Q 



, , dU dU 
~dt~'~dt 



(28) 



Applying first law of thermodynamics in the above equa- 
tion one identifies heat dissipated to the bath as 



dt 



9 



i[/'-^(2a-l)C/'' 



(29) 



Note that, for Stratonovich convention this leads to the 
conventional definition of heat 



Q 



dtxU', 



as expected. Hence from ([25)) and p8| we get 



„/3Q 



P[X\xo] 

PlXlXr] 



(30) 



(31) 



This is the same form as obtained from underdamped 
case and we find that although the Langevin equation 
contains extra terms for the overdamped case in presence 
of multiplicative noise, the modified chain rule ensures 
that the Crooks theorem remains unaffected. Given this 
theorem other fluctuation theorems follow as derived in 
earlier section. 

As in the underdamped case we study the nature of 
probability distribution for work and entropy for simple 
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FIG. 4: Transient work distribution for overdamped case with 
space-independent friction 



model of driven harmonic oscillator for both space in- 
dependent and space dependent case and verifying FTs 
numerically using Heun's method (which is equivalent to 
follow Stratonovich description as discussed earlier) . The 
corresponding Langevin equation is given by 



-y'ix) 
j(x)i = ~kx + Asm(ujt)--^--^T- 

27 



v^^iWm- (32) 



In fig. (Ill) and fig. ([5]), we have plotted the transient work 
distributions for forward and reverse processes, for space- 
independent and space-dependent friction, respectively. 
The functional form of 7(0;) are same as studied in un- 
derdamped case. All units are in dimensionless form and 
taken fc = 1, 7 = 1. We find that the distributions are 
Gaussian for space independent case while for space de- 
pendent it is non-Gaussian. But, the crossing point is 
same. From (e"'^'^) = e"^'^^, the numerically obtained 
free energy difference AF — —0.045, which is equal to 



III. CONCLUSION 
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FIG. 5: Transient work distribution for overdamped case witii 
space dependent friction 
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In this work, we have considered the validity of FTs 
in presence of space-dependent friction, for both under- 
damped and overdamped limit. We find that, although 
no conceptual difficulties arise when the system is under- 
damped, the derivation of the FTs are more involved for 
overdamped system. In latter case, the Langevin equa- 
tion contains extra terms and a careful treatment is re- 
quired to ensure that the Crooks theorem is not con- 
tradicted. This essentially requires invoking the modi- 
fied chain rule [30]. Thus, we conclude that the FTs re- 
main valid for the case of dynamics of a particle in space 
dependent frictional medium, even under coarse grain- 
ing, i.e, reducing the description of the system of two 
phase space variable (x,v,underdamped case), to a sin- 
gle phase space variable (x, overdamped case). Dynam- 
ics of the overdamped case has to be treated with care, 
definition of heat gets modified in different prescriptions 
(Ito, Stratonovich etc) in overdamped case. Proper use 
of generalized chain rule has to be invoked and correct 
propagators for forward and backward processes have to 
be used ^. 

As an illustration, we have analyged the nature of 
P{Astot) and P{W) for the simple case of a driven har- 
monic oscillator in presence of space dependent friction, 
both in underdamped and the overdamped regime (in 
Stratonovich prescription). Distributions for constant 
friction is compared with that of a particle moving in 
a space dependent frictional medium. Moreover, several 
FTs have been verified. This model system is amenable 
to experimental verifiction |31j . 
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In fig.® we have plotted the distribution of total en- 
tropy production for the same example but the dynam- 
ics is overdamped. We found, for space independent 
case the distribution is Gaussian. This is true only for 
equilibrium initial distribution. We have verified sepa- 
rately for initial nonequlibrium distribution P{Astot) is 
non Gaussian. But for space dependent case P(Astot) is 
non-Gaussian even initial equilibrium distribution. Nu- 
merically we find (e^^"'"') = 1.002 which is well within 
our numerical accuracy. 
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